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Abstract 

We explore the transitions B(0 + ) to B tv and B s (0 + ) to B K from lattice 
QCD with Nf — 2 flavours of sea quark, using the static approximation for 
the heavy quark. We evaluate the effective coupling constants, predicting 
a B(0 + ) to B n width of around f60 MeV. Our result for the coupling 
strength adds to the evidence that the B s (0 + ) meson is not predominantly 
a molecular state (BK). 

1 Introduction 

The interest in the excited states of heavy-light mesons has been enhanced by 
the striking discovery that the cs states with J p = + and 1 + have very narrow 
widths [1, 2, 3, 4], This raises the question of whether the corresponding bs 
states will also be narrow. The main reason for the narrow width of D s mesons 
is that the transition to DK is not energetically allowed (for the 2317 MeV 
state) or the state is close to threshold (for the 2457 MeV state) . Thus the only 
allowed hadronic decay proceeds via isospin-violation (since rrid ^ m u ) to D s 7r 
and will have a very small width. Likewise, if the equivalent bs states are close 
to or below the BK threshold, then they will be very narrow. 

Lattice studies have addressed the energies of these P-wave 6s states [5, 6, 7] 
and concluded that they indeed lie close to or below threshold and hence have 
very small decay widths. Although the lattice studies use bs creation operators 
for these states, it is also possible that a molecular description (as a BK bound 
state) is more appropriate, as has been suggested for the cs case [8]. To clarify 
this situation further, it would be very useful to evaluate the hadronic transition 
strength from the scalar B state to a B meson plus a light pseudoscalar meson. 

Here we evaluate these hadronic transition amplitudes using lattice meth- 
ods. This has relevance to the decay of a scalar B or B s meson to B plus a 
pseudoscalar meson. 
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2 Spectrum 



In the heavy quark limit, the Qq meson, which we refer to as a 'B' meson, will be 
the 'hydrogen atom' of QCD. Since the meson is made from non-identical quarks, 
charge conjugation is not a good quantum number. States can be labelled by 
L±, where the coupling of the light quark spin to the orbital angular momentum 
gives j q = L ± 2 ■ In the heavy quark limit these states will be doubly degenerate 
since the heavy quark spin interaction can be neglected, so the S state will have 
J PC = 0~, 1~ while the P_ state will have J p = 0+, 1+. Note that in the static 
case, the self energy of the static heavy quark is unphysical, so that only mass 
differences are physical. Our notation for these static-light mesons is B(nL±) 
where n = 1 (often omitted) is the ground state, and n = 2 the first excited 
states, etc. Here we will be studying the transition from the P_ state to the S 
state emitting a pion in a relative S-wave. This can be applied to the decays 
B(0+) -► B(0-)tt and to B(l+) -► B(l~)ir. 

We shall be using the Nf = 2 lattice configurations [9] with (3 — 5.2 and 
volume 16 3 x 32 with SW-clover improvement coefficient 2.0171. We only use 
the unitary points, namely those with valence light quarks of the same mass as 
the sea quarks. The details of the spectrum from ref. [7] are collected in Table 1. 

The method we shall use to obtain 3-point correlations (next section) using 
timeslice random sources can be used for 2-point correlations and compared 
with the maximal variance reduction (MVR) method [5] used for the 2-point 
correlators in extracting the spectrum [7]. For our lighter quark mass, we find 
the local-local B(S) correlator is more precisely determined for t > 4 by 40 gauge 
configurations of MVR than 100 gauge configurations of time-slice evaluation, 
although the latter had a somewhat smaller computational overhead. Since 
larger t is important for separating ground states and excited states, MVR is 
the method of choice for the 2-point correlation. Because it does not generalise 
efficiently to the 3-point correlation, we use the timeslice method there. 

3 Decay transitions 

Following the methods [10, 12] used to study the hadronic transitions such as 
hybrid decay and p to irir, we can determine the transition amplitude provided 
that there is approximate equality of energies between the initial state B(1P_) 
and the final two-body state B(1S) n. Here we are taking the b quark as static 
and using two flavours of light quark. Staying within the fully unitary sector 
of the theory, we can study transitions with the same valence quarks in the B 
mesons and pion as in the sea. 

The lightest two-body B(1S) it state on a lattice will be when the pion has 
relative momentum zero. The energy differences are then given by aAE = 
0.09, 0.19 for light quarks of mass corresponding approximately to 2m s /3 and 
m s for k = 0.1355 and 0.1350 respectively [7]. Especially for the lighter quark 
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K 

MVR gauges 
t slice gauges 
r /a 

r o m(0 '") 


0.1355 0.1350 
40 20 
100 20 
5.041(40) 4.754(40) 
1.48(3) 1.93(3) 


r m(lS) 
rom(2S) 
rom(lP_) 
r m(2P_) 


3.73(8) 3.68(7) 
5.60(14) 5.61(8) 
4.75(6) 4.71(8) 
7.38(9) 7.1(2) 



Table 1: Lattice parameters and results from ref. [7] for the energies of Qq states 
in units of tq for dynamical fcrmions with Nf = 2. The values of ro/a and the 
qq pseudoscalar meson mass are from ref. [9]. Here we set the scale using ro of 
0.525(25) fm. The heavy-light meson lattice energies contain the static source 
self-energy so that only differences are physical. 




Figure 1: Quark diagram evaluated for the transition from B(P ) to B tt. Here 
B(P_) is the static-light meson that comprises degenerate scalar and axial B 
mesons. 

case, this is an energy difference small enough to apply the method [10, 12], 
namely AEt <5, up to large t- values. 

We need to evaluate correlations between B(1P_) at t = and B(1S) it 
at time t. This involves quark propagators between three space-time points. 
This is illustrated in fig. 1. The heavy quark propagator, however, is trivial to 
evaluate: as a product of gauge links with a (1 ± ji)/2 projector for spin. We 
create B(S) as Q^q and B(P_) as Qq, and, in both cases, also two different 
fuzzed versions of these [5, 7]. In this exploratory study, we only consider a 
pion with zero momentum (so we sum over relative spatial position) with a 
local creation operator qj^q. 

To gain sufficient statistics for the three point correlations, we wish to eval- 
uate the correlation using every space and time point on the lattice as a source. 
To achieve this, we follow the stochastic technique used previously [13, 7]. We 
use a stochastic source ^(complex gaussian random number in every colour, 
dirac, space component) at a given time slice t. We then evaluate the propa- 
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gator (f> from this source using M<j> = £ where M is the Wilson-Dirac matrix 
for the light quark. The required correlation can then be obtained from this 
propagator, schematically as, 

C 3 (t' -t) = <j>*(t',y)ai<i>(t,y)b j (i + i4) j i {Tlt"=t,t>-iU(yX)) ab (i) 

where all repeated indices (and y) are summed. The product of stochastic 
sources implicit in the product of 0's then gives an expectation value which is 
the required sum over pion sources throughout the time-slice at t, whereas the 
noise terms average to zero. By using more independent stochastic samples, the 
average over them (from eq. 1) will have reduced noise. One can also improve 
the signal to noise ratio by combining results from different time values t for 
the stochastic source. 

In practice we used one stochastic sample per time slice, but all time slices in 
turn, as in ref. [12]. This implies 32 inversions per gauge to evaluate the required 
3 point correlation from all sources to all sinks. This is computationally very 
efficient and provides sufficient precision, as shown below. We use 100 gauges 
for the lighter quark mass and 20 for the heavier. 

The motivation for using a source restricted to one time-slice is to ensure that 
the noise contributions decrease as the signal decreases with increasing \t' — t\. 
Note that a two point correlator, for example for the pion with zero momentum 
and local creation and destruction operator, can be obtained likewise from 

C w (t' -t) = <j>*(t',y) ai <j>(t',y) ai (2) 

In this work, we extract the ground state pion contribution to CV from a fit 
to pion correlations obtained from conventional analyses [9] with non-stochastic 
sources, so we do not use the stochastic result for the 2-point pion correlators, 
other than as a check. 

Likewise, for B(P_) with local creation and destruction operators, the two- 
point correlator can be obtained from 

C B( p_)(t' -t) = <Ptf,y) ai Z(t,y) bj (l+> U ) ji (U t „ =t , t ,^U{y,t")) ab (3) 

As we discussed above, this latter expression for C B is more noisy than the 
MVR method at larger \t' — 1\, so we again only use it as a cross check. 

The normalised transition amplitude 2 on a lattice can then be obtained 
from the ratio 

w = = xt + const (4) 

y C B (p_ ) (t)C„ (t)C B (s) (t) 

provided that the transition rate is not too large, namely st < 1. This ratio 
for the decay to n + is plotted for our lighter quark mass in fig. 2. As well as 
illustrating the result for each of our three operators to create a heavy-light 
meson, we can choose to improve the ground state projection of the B(S) and 
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Figure 2: Normalised three particle correlator versus t/a for k=0.1355. The 
points marked (+, x , *) are for local, lightly fuzzed and heavily fuzzed operators 
respectively. The combination which optimises the ground state is shown by 
squares, and a linear fit to it is shown. 



B(P_) by using an appropriate linear combination of local and fuzzed operators. 
The ratio for this improved projection is also illustrated. 

This result shows a linear behaviour, as expected if excited state contribu- 
tions are not significant. We can then read off the hadronic transition amplitude 
ax from the slope - obtaining ax = 0.028(3). 

This is the transition with lattice normalisation and for one quark diagram. 
For the transition B(P_) — > B(S) ir, there will be two quark diagrams contribut- 
ing (since either auord quark pair can be produced, yielding ir + or 7r°) and the 
overall rate will be 3/2 that evaluated from the amplitude x above. To derive 
the appropriate normalisation [10], consider the decay width, even though the 
decay is not energetically allowed with our parameters. Then V = 2irx 2 p where 
p is the two body phase space which evaluates to p = L 3 kE 7r /(2-K 2 ). We have 
an isotropic decay (S-wave) and it is reasonable to assume that x is independent 
of the decay momentum k. 

To increase predictive power, we evaluate the coupling constant in an effec- 
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Figure 3: Normalised three particle correlator versus t/a for k=0.1350. The 
combination which optimises the ground state is shown together with a linear 
fit. We also show the result of a two state model with the correct energy 
difference (aAE = 0.19) between the two body and scalar B states. 



tive lagranian for the three point vertex describing the decay. This coupling 
constant has the dimensions of mass, so we set the scale using the mass of the 
decaying meson. Then the coupling constant squared is proportional to T/k, 
and we use that definition as an effective coupling strength. We then obtain an 
effective coupling strength of 

T/k = 3(L/af(ax) 2 aE v /(2n) = 0.46(9) (5) 

where for the pion at zero monentum, we may use its lattice mass [9, 11] for 
E w . 

In order to explore the dependence on the light quark mass, we use another 
quark mass, although we are limited by the need to keep the transition approx- 
imately on mass-shell, so with decay products of similar energy to the initial 
scalar meson. We used k = 0.1350 where the quark mass is approximately 
strange. We find a similar plot (fig. 3) of xt versus t/a with a slope of 0.0237. 
Since in this case we have a somewhat bigger mismatch (namely aAE = 0.19) 



6 



between the energies for the two body state and of the scalar meson, we can cor- 
rect for this by using a two state model [12]. This shows that we would expect 
some small curvature, even in the ideal case when there is no contribution from 
excited states. Our lattice result is quite consistent with this curvature. Be- 
cause of the additional analysis needed to cope with the larger energy gap, the 
systematic errors are relatively larger in this case. We estimate ax = 0.024(4). 
Then the effective coupling strength is Y/k — 0.46(9), exactly the same value 
as obtained at the lighter quark mass. 

Since a quark- antiquark pair is created in the decay, it might be expected 
that the amplitude to produce heavier quarks was smaller. However, a major 
component of the transition amplitude may come from considerations of the 
overlap of the initial and final states, and this does not depend on the light 
quark mass in any very simple way. Indeed in our study [12] of p decay to two 
pions, we saw some evidence that the decay amplitude was largely independent 
of the light quark mass. This is what we find here for scalar decays also. 

The method we have used to evaluate the hadronic transition is only approx- 
imate, and assumes that the transition amplitude x is relatively small. Since 
we find that xa « 0.03, this is indeed justified. In general, however, one can 
proceed in a rigorous way. This involves determining the energy of the two 
body system (B(15) + n) as accurately as possible with a full QCD lattice 
simulation, and then obtaining the dependence of this on the lattice spatial 
volume. This then gives information on the scattering phase shift in the two 
body channel [14, 15, 16, 17]. It is, of course, consistent to treat the static 
quark as quenched, but all the light quarks need to be treated dynamically. In 
our approach with Nf — 2 flavours of degenerate sea quark, this would allow a 
study of the transition from B(1P_) to B(15) + tt. To have the most accurate 
determination of the two-body energy, one should use a variational approach 
with both two body and one body operators. This will involve the three body 
correlation we have measured above, but also the two to two and box quark 
diagrams. Our preliminary study indicates, as was found for the case of p de- 
cay [12], that these diagrams are too noisy to yield sufficiently accurate results, 
even measuring from all space-time sources for 100 gauge configurations. 

4 Discussion 

As discussed above, we are able to measure the transition amplitude from the + 
bq meson to B tt, provided that the light quark masses are such that the initial 
and final states have very similar energies. For the case we have explored, 
with Nf — 2 flavours of degenerate light quark, this implies that we must 
extrapolate in the light quark mass to make contact with experiment. This we 
do by assuming that the coupling constant for the transition, as described by 
an effective lagrangian, is independent of the light meson mass. This leads to 
the assumption that the effective coupling strength T/k introduced above will 
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be independent of the light quark mass. We do indeed see some evidence from 
our lattice results that this is the case. Thus we shall use our lattice results 
for the reduced width, evaluated where no decays are allowed, to compare with 
experiment and to make predictions. Since we work at a fixed lattice spacing, 
we are unable to estimate the systematic error arising from not taking the 
continuum limit. 

There is a state known experimentally [18] which is a candidate for the + 
bq meson, namely the B** with mass 5698(8) MeV and width 128(18) MeV. 
This corresponds to an effective coupling strength of T/k = 0.34(5). However, 
the experimental state may be a superposition of several states, so mass values 
and widths for the + state are not really known experimentally. 

From lattice studies with static quarks, the excitation energy of the scalar 
B** state is estimated [7] to be 368 ± 31 MeV, where this energy difference 
was evaluated for strange light quarks, but was expected to be similar for non- 
strange light quarks. Using this central value of 368 MeV for the energy release, 
the width of the scalar B** state, with decay to Bw, would be 162(30) MeV. Our 
result is significantly lower than that obtained [19, 20] using a chiral symmetry 
between the ± B mesons, namely a width of around 500 MeV using Ga ~ 1- 

It may also be relevant to compare with experimental data on decays of 
heavy-light mesons with charm quarks, since there is a wider range of data 
available [1, 2, 3, 4, 18]. From the observed [4] mass of 2308±17 ± 15 ± 28 
MeV for D(0+) and width of 276±21 ± 18 ± 60 MeV for decay to D(0") + tt, 
one gets T/k = 0.73^. This is a somewhat larger effective coupling strength 
than the value of 0.46(9) that we obtained above (but consistent within errors), 
although our evaluation is for static quarks whereas charm quarks are known to 
be sufficiently light that this can be a poor approximation for them. 

It is also possible to extract an effective coupling strength for the decay of 
K(1412) to K n, obtaining [18] T/k = 0.48(5). Thus the experimental data are 
consistent with an effective coupling strength of about 0.5 for decays of scalar 
heavy-light mesons with heavy quarks that are b, c and s. This is very consistent 
with our ab initio evaluation which gives around 0.5 also. 

For the bs excited mesons, in the limit of degenerate u and d quarks, there 
will be no decay to pions and the main hadronic decay will be to BK with the 
emission of a light quark-antiquark pair. In this case our evaluation is partially 
quenched, in the sense that the strange quark in the B s meson and K meson is 
not present in the sea. For the decay of a scalar B s meson, the energy release 
may be small or the state may even be stable [7]. Even if the state is stable 
under strong interactions, we can still evaluate the hadronic transition strength 
as an effective coupling. Consider the transition B S (P_) — > B(S) K, there are 
again two quark diagrams, now with equal weight. Our result is then that 
T/k = 0.61(12). If this scalar meson does lie above threshold, we predict a 
width given by that expression. 

Consider now whether the B s (P_ ) meson is a quark-antiquark state or a BK 
meson. Since we have found a non-zero transition amplitude (our x) on a lattice 
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it follows that the meson and the two-body state mix. Indeed when the meson 
mass is degenerate with the two-body energy, there will be an avoided level- 
crossing, with full mixing. What is more significant, however, is the situation 
in a large volume, when the two body energy spectrum becomes continuous. 

The situation in lattice studies is then more like in experiment - one has 
to deduce the composition of a hadron from its observed properties. There 
are lots of extra clues available in lattice studies, however: (i) the mass of 
the state can be explored as the quark mass varies, (ii) the wave-function and 
charge form factor of the state can be measured, (iii) the coupling strength of 
transitions can be evaluated. For the B(P_) meson, lattice studies with Nf = 2 
show a spectrum [7] with a mass which is more-or-less independent of the two 
body (Bit) threshold which would not be expected for a molecular state. The 
coupling strength, as discussed above, is similar to that for the scalar decay 
of K(1412) to K 7r, where a molecular structure is not expected. The charge 
form factor has only been measured [21] for Nf — 2 for the ground state B(S) 
although quenched results [5] show a Bethe Salpeter wavefunction for B(P_) 
similar to quark model expectations. A more definitive lattice conclusion must 
await studies with lighter quarks, but all the evidence at present points to the 
heavy-light scalar meson as not being predominantly a molecular state. 
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